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We study fidelity decay by a uniform semiclassical approach, in the three perturbation regimes, 
namely, the perturbative regime, the Fermi-golden-rule (FGR) regime, and the Lyapunov regime. A 
semiclassical expression is derived for fidelity of initial Gaussian wave packets with width of the order 
\/fL {h being the effective Planck constant). Short time decay of fidelity of initial Gaussian wave 
packets is also studied, with respect to two time scales introduced in the semiclassical approach. 
In the perturbative regime, it is confirmed numerically that fidelity has the FGR decay before the 
Gaussian decay sets in. An explanation is suggested to a non-FGR decay in the FGR regime, which 
has been observed in a system with weak chaos in the classical limit, by using the Levy distribution 
as an approximation for the distribution of action difference. In the Lyapunov regime, it is shown 
that the average of the logarithm of fidelity may have roughly the Lyapunov decay within some 
time interval, in systems possessing large fluctuation of the finite-time Lyapunov exponent in the 
classical limit. 
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I. INTRODUCTION 

It is well known that, in classically chaotic systems, 
time evolution of trajectories in phase space is sensitive 
to small changes in initial conditions, as well as in pa- 
rameters in the Hamiltonians. On the other hand, time 
evolution of state-vectors in Hilbert space is insensitive 
to small change of initial conditions. Nearly twenty years 
ago, Peres observed that small change in perturbation 
parameters can be employed in the study of the stability 
of quantum motion in the Hilbert space 0, supported 
by further numerical investigations 0, 1^ ■ The quantity 
used to measure the stability of quantum motion is the 
quantum Loschmidt echo, or fidelity in the field of quan- 
tum information ^|^. It is an overlap of the evolution of 
the same initial state under two Hamiltonians with slight 
difference in the classical limit, M{t) — \m{t)\'^ , where 

m{t) = («'o|exp(ii7ii/?i)exp(-ii7oV^)l*o)- (1) 

Here Hq is the Hamiltonian of a classically chaotic system 
and Hi — Hq + eV with e being a small quantity. 

Fidelity decay has attracted increasing attention, since 
the work of Jalabert and Pastawski 6], which relates 
the decay rate of fidelity to the (maximum) Lyapunov 
exponent of the underlying classical dynamics. In or- 
der to understand the behavior of fidelity in various 
systems , ex tensive investi gati ons have been carried out 

lilBSlillillilM 

illMillMllllMlllllg- Previous investigations 
show the existence of at least four regimes of perturba- 
tion strength for fidelity decay : (i) In the pertur- 
bative regime below a perturbative border, fidelity has a 
Gaussian decay. In this regime, the typical transition ma- 
trix element is smaller than the mean level spacing, (ii) 
Above the perturbative regime is the Fermi-golden-rule 
(FGR) regime, with an exponential decay of the fidelity, 
M{t) (X exp(— Pi), where P is the half-width of the local 



spectral density of states (LDOS) M- The decay rate can 
also be calculated semiclassically [13 • (iii) With increas- 
ing perturbation strength, one enters into the Lyapunov 
regime, in which M(t) cx exp(— At), with A being the Lya- 
punov exponent of the underlying classical dynamics jj] . 
(iv) In the regime above the Lyapunov regime, perturba- 
tion strength is so large that the classical perturbation 
theory fails. Presently, little is known about the decay- 
ing behavior of fidelity in this regime. Within a random 
matrix theory approach, in which fidelity is expressed as 
the Fourier transform of LDOS, a Gaussian decay was 
suggested for fidelity decay in this regime , without 
further numerical confirmation (cf. |3lL l32| for shape of 
LDOS in this regime). It is also known that, for time t 
short enough, fidelity has a quadratic decay, which may 
be extended to a Gaussian decay, just as a direct result 
of perturbation theory . 

Most recent investigations show that the above picture 
of fidelity decay is incomplete, at least in four aspects. 
Firstly, in the perturbative regime, numerical results p3 | 
show that fidelity in the kicked rotator model has an 
exponential decay, which can be described by their semi- 
classical approximation, before the Gaussian decay sets 
in at about the Heisenberg time. A random matrix ap- 
proach to fidelity also suggests an approximately expo- 
nential decay of fidelity at t short enough ^.29] • It is not 
quite clear whether this exponential decay is the FGR 
decay or not. 

Secondly, a non-FGR decay of fidelity in the expected 
FGR regime has been found in a system with weak chaos 
in the classical limit, which is induced by deviation of 
the distribution of action difference from the expected 
Gaussian distribution ^28J. An analytical description for 
the rate of the non-FGR decay is still lacking. 

Thirdly, in the Lyapunov regime, the decay rate of av- 
erage fidelity has been found different from the Lyapunov 
exponent, although still perturbation-independent, in 
systems possessing large fluctuation in the finite-time 
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Lyapunov exponent |33. l3J[. as in the kicked top and 
kicked rotator models Il9l |20| . A semiclassical WKB 
description of wave packets suggests an exp(— Aii) decay 
for the fidelity, with Ai < A |23|. More recently, more 
general semiclassical expressions of fidelity decay have 
been derived, with the Lyapunov decay and the Ai decay 
being two limiting cases oU\ , along the lines of the semi- 
classical treatment to fidelity in Refs. 0, IM lIOj 22, 28] . 
However, the situation has not been clarified completely, 
since numerical results in the kicked top model show that 
the Lyapunov decay can be resumed in an approximate 
way, if average is performed on the logarithm of fidelity, 
but not on the fidelity itself (191] . 

Finally, in the deep Lyapunov regime, fidelity of initial 
Gaussian wave packets may have a decay which is super- 
exponential and much faster than the Lyapunov decay 
at short initial times Meanwhile, a decay with a 

rate of twice the Lyapunov exponent may appear before 
a time scale introduced in [sil ] , in systems with constant 
local Lyapunov exponents. A quantitative description of 
the former decay is still not available and the time scale 
that separates the two faster than Lyapunov decays is 
unknown 

In this paper, we use the uniform semiclassical ap- 
proach introduced in |22l ] to study the problems men- 
tioned above. This approach is not only a suitable 
method for numerical evaluation of fideli ty, but also a 
good starting point for analytical study jl^, EO]- For 
simplicity, we study ID kicked systems only in this pa- 
per. 

The paper is organized as the following. In Scct.m we 
introduce two models, the kicked rotator and the saw- 
tooth map, that will be employed for numerical check 
of our analytical results. A major difference between 
the two models is that the sawtooth map has a constant 
finite-time Lyapunov exponent, while the kicked rotator 
has large fluctuation in the finite-time Lyapunov expo- 
nent. For the sawtooth map, the semiclassical prediction 
of the rate of FGR-decay can be calculated accurately 
at some parameter values, meanwhile, it has weak chaos 
with a structure of cantori in some parameter regime |35j] . 

The validity of the uniform semiclassical approach has 
been checked numerically for initial point sources p3 . l2^ . 
For initial Gaussian wave packets, narrowness of the 
packets are assumed in deriving the semiclassical expres- 
sions of fidelity in "g*, "2^. In Sect. IIIll we show that the 
expression in J2,] fails in describing the fidelity decay, 
when the width of the initial Gaussian wave packet is of 
the order Vh, where h is the effective Planck constant in 
the ID kicked systems studied here. By considering the 
second order term in the Taylor expansion of the action, 
we derive a modified expression, which works well for this 
kind of initial wave packets. 

In Sect. II Vl we discuss short time behavior of fidelity. 
In particular, a time scale is introduced for fidelity de- 
cay of initial Gaussian wave packets, which separates the 
two faster than Lyapunov decays mentioned above, and 
an analytical expression is derived for fidelity before this 



time scale. Dependence of the first-kick decay of fidelity 
of initial point sources on perturbation strength is also 
derived. Fidelity decay in the perturbative and the FGR 
regimes are studied in Sect. with emphasis on the 
problems mentioned above. In the perturbative regime, 
the exponential fidelity decay before the Heisenberg time 
is shown numerically to coincide with the FGR decay 
in the sawtooth map. In the FGR regime, we revisit 
the non-FGR decay in the sawtooth map found in [28| . 
and show numerically that the central part of the distri- 
bution of action difference can be approximated by the 
Levy distribution, which can explain some properties of 
the non-FGR decay. Section IVTI is devoted to a study of 
fidelity decay in the Lyapunov regime, when the average 
is performed over the logarithm of fidelity. Conclusions 
and discussions are given in Sect. lVIll 

II. MODELS: KICKED ROTATOR AND 
SAWTOOTH MAP 

The Hamiltonians of the two models employed in this 
paper are of the forms, 

^ oo 

^ - 2^' + ^^^'^ ^""^ ^ ~ ""^^ ^'^^ 

Vfe(r) — K cosr for kicked rotator, and (3) 

Vs{r) = -K{r - iTf/2 (for sawtooth map). (4) 

For simplicity, the period T is set to be unit, T ~ 1. Kicks 
are switched on at t = n, n = 0,1,2,.... The classical 
map describing the kicked rotator is the standard map, 

Pn+i ^Pn + K sin(r„) (mod 27r), 
rn+i=rn+Pn+i (mod27r). (5) 

The sawtooth map is 

Pn+i = Pn + K{rn- tt) (mod 27r) 
rn+i=rn+Pn+i (mod27r). (6) 

Equation ((BJ can be rewritten in the matrix form. 




where the 2x2 constant matrix possesses two eigenval- 
ues 1 + {K ± +4X)/2. At K > 0, motion in the 
sawtooth map is completely chaotic, with the Lyapunov 
exponent \ = ln{{2 + K + [(2 + Kf - 4]i/2)/2}, given 
by the largest eigenvalue of the matrix. The finite-time 
Lyapunov exponent has the same value as the Lyapunov 
exponent A defined in the limit t ^ oo. On the other 
hand, the standard map, which is chaotic at K larger 
than 6 or so, does not have a constant finite-time Lya- 
punov exponent, because the mapping matrix is a func- 
tion of r„. It is of interest to mention a recent result on 
the classical analog of fidelity, namely, for systems with 
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more than one dimensional configuration space the clas- 
sical fidelity has a decaying rate related to not only the 
maximum Lyapunov exponent, but also other positive 
Lyapunov exponents |36| . 

The two classical systems are quantized on a torus 
[37! l3^ l3^ 23 ■ ^ system with ID finite configura- 
tion space, < r < Tm, and ID finite momentum space, 
< p < Pm, the effective Planck constant h and the 
dimension N of the Hilbert space has the relation 



Nh = rmPr) 



(8) 



In both models, we take = p,„ = 27r, hence, T% = 
2t:/N. 

Floquet operators in the two quantized systems have 
the form 



U = exp[-ipV(2?i)] eM-iVk(s){f)/K\- 



(9) 



Eigenstates of f are denoted by \ and f\j) ~ j), with 
j = 0, 1, 2, . . . , iV — 1. In this representation, elements of 
the operator U are 



1 



N 



exp 



N 



27r 



(10) 

The evolution of states, = U^ipo, is calculated nu- 
merically by the fast Fourier transform (FFT) method. 

The fidelity in Eq. ^ involves two slightly different 
Hamiltonians, Hq and Hi = Hq + eV. In what follows, 
Hq takes the form of H in Eq. (|2Jl , and 



except in Sect. lIVB2l 
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(11) 
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III. UNIFORM SEMICLASSICAL APPROACH 
TO FIDELITY 

A. Approximation to fidelity witli action expanded 
to the first order term 

For the sake of completeness and convenience in pre- 
senting our results, we briefiy recall the main results of 
the semiclassical approach to fidelity in Refs. 0, in 
this section. 

In the semiclassical approach, an initial state ipoi'^^o) in 
a d-dimensional configuration space is propagated by the 
semiclassical Van Vleck-Gutz wilier propagator, 



(12) 



'0sc(r;t)= J droA'sc(r,ro;i)V'o(i'o), 
where Ksc{r,ro;t) = Y.s ^^(^^^o;t), with 



a 



1/2 



rexp 



1 17T 



(13) 



Here, the label s [more exactly s(r,ro;t)] indicates clas- 
sical trajectories starting from Tq and ending at r in a 
time t, the action Ss{r,ro;t) is the time integral of the 
Lagrangian along the trajectory s, Ss{r,ro;t) = j^dt'C, 
and Cs = \det{d'^ Ss/droidrj)\. fig is the Maslov index 
counting the conjugate points. 

Consider an initial Gaussian wave packet centered at 
To, with dispersion ^ and mean momentum po. 



d/i 

^o(ro) = ( ) exp 



1 



(14) 



When 5 is small enough, within the effective domain of 
ro, Ss{v,ro;t) can be expanded in the Taylor expansion 
with respect to the center Tq, 



Ss{r, ro; t) = Ss{r, ro; t) - (ro - ro) • 
where 

dSs{r,ro;t) 



Ps = 



(15) 



(16) 



is the initial momentum of the trajectory s(r,ro;i)- 

The semiclassical approximation to the fidelity ampli- 
tude m{t) in Eq. (P) is 

mit) ^ Jdr [V4^ (r; t)] * i^S:" (r; t), (17) 

where the two states are propagated by the two Hamil- 
tonians Hi and Hq, respectively. For quite small ^, the 
expansion in Eq. IjlSI) can be truncated at the first order 
term. Then, by using Eqs. (|13|I - H16() the integration on 
the right hand side of Eq. H12|l is calculated. 
The amplitude m{t) thus obtained is (|, 



m{t) 
dr CsCxp 



msci(i) = ( 

nh 



d/2 



-ASs{r,ro;t)- ^{p, 
n n 



Po)' 



,(18) 



where A5s(r,ro;t) is the action difference for the two 
trajectories with the same label s in the two systems Hi 
and Hq. In the first order classical perturbation theory, 
the difference between the two trajectories s is assumed 
negligible, 

ASs{r,ro;t) ^ e f dt'V[rit')] (19) 
Jo 

with V evaluated along the trajectory. 

A simpler expression of m(t) can be obtained by chang- 
ing the variable r — > po 22] , 



nisciit) = ( 



d/2 



dpoexp 



* A ~ ,x (PO - Po)^ 

^A5(po,ro;0-^^y^ 



(20) 
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where AS{po,ro;t) coincides with AS's(r,ro;i) for the 
same trajectory s with initial momentum pp. The main 
contribution to the right hand side of Eq. Ij20|l comes 
from a window in the po space, which is centered at po 
and has a size of the order h/£^. 

For a system with finite momentum space, Eq. H20|l is 
invahd for initial Gaussian wave packets that are wide 
in the momentum space. The extreme case is for initial 
point sources, {r\"^o) = \/{'2nh)'^/Vp6{r — Tq), with Vp 
being the volume of the momentum space. In this case 



m{t) ~ TOp(ro, t) 



1 



dpoexp 



-AS'(po,ro;i) 



(21) 

The above semiclassical expressions of m{t) suggests in- 
troducing 



a = e/?i, 

as a quantum perturbation parameter. 



(22) 



B. Contribution to fidelity from the second order 
term in action expansion 

Hereafter we restrict our discussions to 1-dimensional 
kicked systems. In deriving Eq. H20() for fidelity of an ini- 
tial Gaussian wave packet, the right hand side of Eq. H15|) 
is truncated at the first order term. Hence, Eq. (|20|l is 
valid only when S^- <^h, ov 1 with the parameter k 
defined by 



(23) 



which has been confirmed in our numerical calculation. 

When the condition k 3> 1 is not satisfied, high or- 
der terms on the right hand side of Eq. (|15|1 may give 
considerable contribution. Indeed, numerically, obvious 
deviation of Msci(t) = |msci(i)P from the exact M(t) 
has been observed at n close to one or smaller (see Fig. ^ 
for an example) . We remark that numerical evaluation of 
the right hand side of Ea. (|20|) for msci(i) becomes more 
and more difficult with increasing t, because the num- 
ber of the oscillation of AS'(po,n];0 versus po increases 
exponentially ^30,] . 

To have a good semiclassical approximation at k 1, 
one needs to consider the second order term on the right 
hand side of Eq. ((T3|l . 

Ss{r,ro;t) ~ Ss{r,ro;t) - (ro - ro)ps - T^^i^'o ~ 

2 dro 

(24) 
• (25) 



where 



dps ^ 9ps 
dro ~ dro 



d^S,ir,ro;t) 



5rg 




FIG. 1: Comparison of the exact values of a single fi- 
delity M{t) and their semiclassical approximations Maci(t) in 
Eq. 11201 and Msc2{t) in Eq. 12611 . in the kicked rotator model. 
Parameters are K ^IQ,N = 2^^ = 131072, cr = 1, and k = 1. 
The typical value of D in Eq. 1271 is about 8.9 ^ 1, ex- 
plaining the deviation of Msci (t) from the exact values, while 
Msc2(i) ~ |w-sc2(t)P being quite close to the exact ones. 



Using Eq. (|24|l and following a procedure similar to the 
derivation of Eq. H2UII , we obtain 



dpo 



msc2it) 



0F no 



exp 



-AS{po,ro;t)^^^^ 



,(26) 



where 



D = 



(27) 



Note that D is a function oipo, ro and t. A numerical test 
for this modified semiclassical approximation is shown in 

Fig.n 

Comparing the two equations H2Q|I and H26(l . it is seen 
that the modification is to replace {h/^) in Eq. (pn|l by 



Wr, 



(28) 



i.e., the change in the size of the effective window for inte- 
gration. Therefore, the modified semiclassical expression 
in Eq. (|26|l predicts the same long-time decaying behav- 
iors of fidelity as Eq. (|2()|l , more precisely, the same decay- 
ing rate for the FGR decay in the FGR regime, and the 
same Ai{t) decay in the Lyapunov regime (cf. Sect. lVip . 

When the value of k decreases further, higher order 
terms in the Taylor expansion of the action should be 
considered, e.g., at ~ Ti, the third order term should 
be taken into account. 
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FIG. 2: Variation of AS{po,ro;t) /2-Ke versus po/2-K in the 
sawtooth map at K = 1, for a value of ro taken randomly 
within [0,27r). r(t) is the position r at time t, with initial 
condition (ro,po)- For clarity, r is plotted as a continuous 
function of po, by adding 2nn at the discontinuous points. 



IV. SHORT TIME BEHAVIOR OF FIDELITY 

The behavior of fidelity in short times is initial-state- 
dependent. In this section, we discuss short-time-decay 
of fidelity of initial Gaussian wave packets and of initial 
point sources. 



A. Oscillation of AS{po,ro;t) versus po 

The semiclassical expressions of fidelity discussed in 
the previous section, specifically, Eqs. (^0)) . and 
show that decaying behavior of M{t) is mainly deter- 
mined by the action difference AS{po, rp; t) as a function 
of pq. Therefore, before addressing fidelity decay, it is 
useful to first discuss an important property of the ac- 
tion difference, namely, its oscillation as a function of pq. 

The number of the oscillations of AS", as po runs over 
[0, 27r), increases exponentially with time t. To see this, 
using Eq. H19|) . we write the slope of AS'/e, denoted by 
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FIG. 3; Same as Fig.|5|but for the standard map of K = 10. 



fcp, in the following explicit form. 



l dAS{po,ro;t) 
e dpo 



' ^^,dV dr'{t') 
dr' dpo 



(29) 



where the dependence of kp on po,ro and t is not writ- 
ten explicitly, for brevity. Due to the underlying chaotic 
classical dynamics, \dr' {t') / dpo\ increases exponentially 
with t' , on average. On the other hand, the variance of 
AS* increases as t hence, the typical value of \AS\ 
increases as ^/i. Therefore, the number of the oscillations 
of AS must increase exponentially with t. The fast os- 
cillation of AS* is crucial for understanding the long time 
decay of fidelity in the Lyapunov regime [2^ |33| • 

We present some examples of the oscillating behavior 
of AS* in the sawtooth map (Fig. [21, as well as some in 
the kicked rotator model (Fig.lSJ. The two models have 
different dependence of the position r{t) on the initial mo- 
mentum Pq at fixed time t. Specifically, in the sawtooth 
map, r is a linear function of po except at the discontin- 
uous points, with the slope given by the constant local 
Lyapunov exponent, while in the kicked rotator, it is an 
oscillating fimction. 



B. Time scales ri and t-2 for fidelity of initial 
Gaussian wave packets 

Fidelity of initial narrow Gaussian wave packets has 
rich behavior at short times. For example, there are both 
very fast and quite slow decays at the first several kicks 
in the deep Lyapunov regime j20j , as well as a decay with 
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a rate of twice the Lyapunov exponent '3Cr| . By using the 
uniform semiclassical approach discussed above, we give 
a unified description for these phenomena in this section. 



1. Time interval t < ti 

The main contribution to the right hand side of 
Eq. H26II comes from a narrow window in the pq space. 
For time t short enough, hnear approximation can be 
used for the action difference AS within the narrow win- 
dow. This suggests the introduction of a time scale, de- 
noted by Tl, such that for t < ti hnear approximation 
of AS* can be used in calculating the right hand side of 
Eq. p)|. 



that A(ti) ~ A, the dependence of ri on h is given by 
(l/2A)lnn"^ for ^ = H^^^, which is half the Ehrenfest 
time. 

When the change of D is negligible within the effective 
narrow window of po, in which the linear approximation 
for AS" in Eq. (pn|l can be used, substituting Eq. (p!n|l into 
Eq. if^ . we have 



Msc2{t) ~ exp 



t<Tu (35) 



with the time dependence on the right hand side given 
by kp(ro,po;t). 

Due to exponential divergence of neighboring trajec- 
tories in phase space, the main contribution to the right 
_ hand side of Eq. (|29|l comes from times t m t' . Therefore, 

AS{po, ro; t) ~ AS'(po, n); t) + ekp{po - pq), (30) \kp\ increases typically as Cfce'^^*)*, with Cfe being the pre- 



where kp is the value of kp in Eq. H29(l at the center 
(ro,_Po) of the initial Gaussian packet. 

To give an estimation to ti, we use Apo(i) to denote 
the size of the region in the pa space, which is capable 
of the above linear approximation for AS*. One should 
note that Apolt) shrinks exponentially, due to the expo- 
nentially increasing of the number of oscillations of AS' 
versus po- Since the oscillation of AS is mainly induced 
by local instability of trajectories, Apo^t) shrinks roughly 
as e"'^'-*-'*, where 



Ait)= _ lim - [ln\6x{t)/Sxm] 

dx{0)—>-0 t 



(31) 



with 6x{t) denoting distance in phase space and aver- 
age performed over phase space. (In a classical system 
with strong chaos, A{t) usually approaches the Lyapunov 
exponent A quickly, as will be illustrated numerically in 
Sect.EIl) Then, 



Apo(i) ^ fo(i)Apo(l)e~^^*^-(*~'\ 



(32) 



where b{t) is the influence of other factors, such as the 
variance of AS" increasing linearly with t, and changes 
much slower than the exponential term. At time ti, we 
write 



Apo{Ti) = aiWp, 



(33) 



where ai > 1 is determined by the accuracy required. 

Substituting Eq. into Eq. for i = n, we ob- 
tain 



1 



A(ri: 



In 



6Apo(l) 



aiWp 



(34) 



where b is the average value of b{t) for small t. Several 
points can be seen in Eq. Firstly, since Apo(l) 

decreases with increasing a, when a is large enough, 
the right hand side of Eq. (|34|) can be smaller than 1, 
implying ti = 0. Secondly, for large enough ti such 



factor. For this time dependence of j/cpj, Eq. H35|l predicts 



Msc2(t) ^ exp 



t<Ti. (36) 



This gives the extremely fast, double-exponential-rate 
decay of fidelity predicted in Ref. M{t) ~ 

exp(— const X e^'^*), when A{t) ^ A. Note that our predic- 
tion H35|l is more general, since it works for small values 
of \kp\ as well, so long as po is not quite close to any 
stationary point, which may invalid the approximation 
in Eq. 

Numerical check of our prediction H35|) is shown in 
Fig. 21 for the first kick. For a < 40, the analytical 
results have good agreement with the exact numerical 
calculations. With increasing a, deviation enlarges, with 
Msc{t) < M{t), because the difference between the exact 
values of the phase AS/h and their linear approximations 
increases linearly with a. 



2. Time interval ti < t < T2 

For t > Tl, the main contribution to the right hand side 
of Eq. (|26|l comes from the integration over the region 
[pq — Wp,po Wp]. It is useful to introduce a second time 
scale, denoted as T2, at which AS{po,ro;t) completes one 
full oscillation period as po runs over [po — Wp,po + Wp] 
|30| . Note that ti < T2, according to their definitions. 

In order to estimate the time T2, we note that the num- 
ber of oscillations of AS* increases as e"^*^*)*, then, T2 sat- 
isfies the following relation 



1 



T2 



■In 



/ TT \ 



When A(t2) 



A(r2) \coWpJ ' 
A, this relation gives the estimation 
1 



T2 



A^N 

A Vco% 



(37) 



(38) 
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FIG. 4: Comparison of the exact values of M{t) and the 
semiclassical prediction given by Eq. 1351 at the first kick 
t = 1, for one initial Gaussian wave packet, in the kicked 
rotator model with K = IQ, n = 1. 



The time scale T2 is important in understanding short- 
time decay of fidelity in the deep Lyapunov regime with 
cr ^ 1. Indeed, in the time interval ri < t < T2, the 
phase {AS/h) on the right hand side of Eq. if^ . as a 
function of po, can usually be approximated by a straight 
line within the region po € [Po — Wp,po + Wp]. Then, for 
initial states satisfying \akp\ '3> n/wp, one has [30l| . 



Tl <t < T2. 



(39) 



To be more specific, let us consider a special kind of 
system, which has constant local Lyapunov exponent A 
in the classical limit and has no stationary point of AS*, 
i.e., kp ^ for all pq. For such systems, A{t) — A. As 
shown in 's^, when the smallest |fcp| are so large that 
Eq. (|39|) is applicable for all initial states, the average 
fidelity has a decay with a rate of twice the Lyapunov 
exponent. 



M{t) 



~2Xt 



(40) 



since kp increases as e 



^* on average. Here we are inter- 
ested in the general situation, in which the smallest \kp\ 
are not large enough for the applications of Eq. (|39|l . i.e., 
when \<7kp\ 3> Tr/wp is not satisfied for \kp\ close to its 
smallest value. In this general case, the average fidelity 
has a decay rate smaller than 2A due to the influence the 
small \kp\. On the other hand, we note that the size of 
the region of po with quite small \kp\ decreases exponen- 
tially, due to the exponential increment of the number of 
oscillation of AS', therefore, the decay rate of the average 
fidelity should be larger than A. 

To verify numerically the above prediction of fidelity 
decay between e""^* and e"^"^* for ti < t < T2, we use the 
sawtooth map, with the following form of the perturba- 



FIG. 5: Same as Fig. |21but for the sawtooth map with per- 
turbations in Eq. I|41|l and t — 2. 



tion V {Ho unchanged) 

oo 

yii) = Yi^) m _ „y) ^here (41) 

ri=0 

V^'\r) = -M^ir - Try , i = 1,2,3,4,5. (42) 

Setting the coefficient A/'2 — 1/2, V^^'' gives the pertur- 
bation in Eq. (|ll|l . which are also used in 0,|2a|- The 
other coefficients Afi are chosen by the requirement of 
having the same decaying rate in the FGR regime, i.e., 
possessing the same value of the classical action diffusion 
constant K{E) QJ^. For kicked maps, K{E) has the form 

MM, 



oo 

K{E) = -C(0) + C{1), where 



(43) 



1=1 



C{1) = {{V[r{l)] - {V)}{V[rm - {V)}) (44) 

with the average performed over phase space. 

At integer values of K in the sawtooth map in Eq. 10}, 
simple derivation shows that C{1) = Q for I ^0 and 



1 



C(0): 

Then, 



2i + l 



(2i + l)(i + l) 



(odd i); 
(even i). 



(45) 



15 



v/1.4 
"3^ 



, Ma 



47r2^ 



37r3 



(46) 



For the sawtooth map with the above perturbation, 
the action difference at time t can be written as 



AS - -e^M[r(n)-7r]* 



(47) 



n=0 



8 




FIG. 6: Values of kp of a rg chosen randomly, in the sawtooth 
map with perturbation V''^ of i = 3 and 5, t = 3. 



It is easy to prove, by using Eq. that r{n) for any 
fixed n is a monotonically increasing function of po, ex- 
cept at the discontinuous changes from to 27r (or re- 
versely). Then, Eq. H47|l shows that no point exists at 



which kp = for odd 



i, while kp can be zero for even i. 



(See Figs. 12 and El for some examples of numerical illus- 
trations.) Therefore, in the time interval ti < t < T2, 
M(t) should have the standard Lyapunov decay for even 
i |3Clj | and faster than Lyapunov decay for odd i. 

Some values of \kp\ for V^^^ of j = 3 and 5 are pre- 
sented in Fig. It is seen that some kp for i = 5 are 
quite close to zero, implying a decay rate of fidelity be- 
tween A and 2A in the time interval (ti,T2) for i = 5; on 
the other hand, the smallest \kp\ for i = 3 is not quite 
close to zero, implying a decay rate of 2A. Indeed, these 
predictions have been confirmed in a direct calculation 
of M{t), as shown in Fig. [T] The values of ti and T2 can 
be estimated as follows. At if = 1 and ^ — Vh, numeri- 
cal computation shows that Cq — 0.45 and D ~ 1.9. We 
take Apo(l) ~ 27r/100 for a = 100 (cf. Fig. |3 for vari- 
ation of AS/e at i = 1), ai 5, 5 1, then, Eq. 
gives Ti ^ 1. Meanwhile, Eq. H38|l gives T2 « 6.5. The 
two estimations are in good agreement with the direct 
numerical results shown in Fig. [71 



C. Dependence of first-kick-decay of fidelity on 
perturbation strength for initial point sources 

Fidelity of initial point sources, described by Eq. 12111 
with integration performed over the whole pq domain, has 
a short-time decay different from that of initial Gaussian 
wave packets discussed above. 

In the FGR regime, fidelity can be calculated by writ- 
ing the right hand side of E g. 12111 in terms of the distri- 
bution of AS" (see Refs. [H l2a|l When the distribution 
of AS is close to the expected Gaussian distribution, one 




2 4 6 8 10 1214 2 4 6 8 1012 1416 
t t 



FIG. 7: Decay of the averaged fidelity in the sawtooth map 
with K =1, for V'-''> of i = 1, 3, 4, 5 in Eq. igTJ, with param- 
eters a = 100, A'' = 131072, ^ — Vh. For these parameters, 
ri ~ 1 and T2 ~ 6.5. In the time interval ri < t < T2, the 
average fidelity has the Lyapunov decay for i = 4, a double- 
Lyapunov decay for i = 3, and a decay between the two de- 
cays for 1 = 5. The quite fast decay for i = 1 is due to the 
linear dependence of AS on po, as shown in Fig. |3 Aver- 
ages are performed over 2000 initial Gaussian packets, with 
centers taken randomly with fiat distribution in the region 
n/2 < fo(po) < 37r/2. 



has the FGR decay for fidelity. 



(48) 



where K{E) is the classical action diffusion constant in 
Eq. 

In the deep Lyapunov regime, Eq. 121(1 enables an es- 
timation to the dependence of fidelity on perturbation 
strength, which is l/cr, as shown in Eq. l|A4p in Appendix 
1X1 Since this dependence does not change with time, it 
can be seen at the first kick only. For systems with con- 
stant local Lyapunov exponents, combining Eq. ((A4p and 
the known Lyapunov decay, we have 



Mp{t) cx - exp(-At). 



(49) 



Figure IHl presents an example of numerical confirmation 
to the above predictions for the first-kick decay of fidelity. 
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FIG. 8: Decay of fidelity M{t) at the first kick t = 1 as a 
function of a, for a single initial point source, in the logarithm 
scale. It shows a FGR behavior for small a and a l/cr depen- 
dence for large a. MYGR{t) — e"^'^"^ ' in the sawtooth map 
with Ko = 1. 
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FIG. 9: Decay of the averaged fidelity in the perturbative 
regime of the sawtooth map ai K = 1, N — 512, and a = 
10~^. Average is performed over 100 initial point sources 
taken randomly in the configuration space. The decay follows 
the FGR decay approximately up to t ~ 400. 



V. THE PERTURBATIVE AND THE FGR 
REGIMES 

In this Section, we study fidelity decay before the 
Heisenberg time in the perturbative regime, and the influ- 
ence of weak chaos on fidelity decay in the FGR regime. 



A. The perturbative regime 

The regime of quite small e, more precisely, quite small 
(T, is named the perturbative (PT) regime, in which fi- 
delity has a Gaussian type decay 0, 0, . Combin- 
ing the perturbation theory, the random matrix theory 
(RMT), and the semiclassical theory, it has been found 
that 

MMt)^e^p(-^^§^ahA, (50) 

for quantized maps, where 2g/ (3 is the number of classical 
orbits with identical action, which is 2 for the models 
used here, and d is the total mean density of states. The 
index f3 — 1 for time-reversal-invariant systems and /? — 
2 for time-reversal-breaking systems. (Cf. Ref. 0| for 
the expression for continuous variables.) 

The Gaussian decay H50|l sets in at about the Heisen- 
berg time tn ~ N. For t short enough compared with 
tn, there is no reason for the semiclassical approach to 
fail in describing the fidelity decay, even for quite small 
a. Therefore, it is reasonable to expect the FGR decay 
in Eq. H48|) for fidelity in the perturbative regime, before 
some time shorter than tn- In fact, numerically, fidelity 
has been found to possess an exponential decay in the 



kicked rotator model, before some time shorter than the 
Heisenberg time Ih 22]. Whether this exponential de- 
cay is the FGR decay or not is unclear, since the value 
of K{E) in this model is calculated only approximately. 

To make the situation clearer, we employ the sawtooth 
map, in which the value of K{E) can be computed an- 
alytically for integer values of the parameter K in its 
Hamiltonian, K{E) = ir^/m ~ 1.08 [see Eq. (gSJ]. Nu- 
merical results indeed support the above argument that 
fidelity has the FGR decay, before a time shorter than 
the Heisenberg time (see Fig. IS)). 

Deviation of fidelity decay from the FGR decay in 
the perturbative regime, provides a good opportunity for 
a numerical study of the breakdown time of the semi- 
classical approach, denoted hy Ib- It is known that 
t B is proportional to some inverse algebraic power of Ti 
|M IM EE m ' In the sawtooth map, it was found that 
Ib is linear in \/fi, more exactly, ts ~ O.SA^ oc Ti~^ , as 
seen in Fig. ^1 

B. The FGR regime 

When the perturbation parameter u is large enough, 
M{t) becomes close to its saturation value, which is pro- 
portional to 1/7V before the Gaussian decay sets in 
at about the Heisenberg time, then, one enters into the 
FGR regime. The perturbation border for the crossover 
from the perturbative regime to the FGR regime can be 
estimated as ^3 



''P-h^^2K{E)N- ^^^^ 
For systems possessing strong chaos in the classical 
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FIG. 10: The breakdown time ts of the semiclassical ap- 
proach versus the dimension A'^ of the Hilbert space, in the 
logarithmic scale for the sawtooth map with K — 1, ai differ- 
ent values of (T. In numerical calculation, is taken as the 
first kick at which the relative error | [M{t) - MFGR{t)]/M{t)\ 
is larger than 0.1, where MpGR(t) is the FGR decay in Eq. 14811 
predicted by the semiclassical theory. 



FIG. 11: Decay rate 7 of fidelity versus perturbation strength 
a, calculated by the best fit of —yt to In M{t), for some values 
of the parameter K in the sawtooth map. In the calculation of 
average fidelity, 1000 randomly chosen Gaussian wave packets 
are used as initial states. The solid curve shows the rate 
r ~ 2.2(7^ of the FGR-decay. iV = 131072, (, ^ Vh. 



limit, two analytical approaches are available to obtain 
the FGR decay, namely, the semiclassical approach with 
the assumption of a Gaussian type distribution of A5'|0 , 
which gives the result in Eq. H48() . and the RMT ap- 
proach making use of the closeness between the form of 
the LDOS and the Lorentzian distribution 0, 01 ■ The 
two approaches are believed to be equivalent, when both 
are valid; while an analytical proof of the equivalence is 
available only in some special cases 0. An interesting 
phenomenon is that the two approaches are complemen- 
tary in some cases. For example, for quite short time 
there is no analytical reason for the distribution of AS' to 
be close to a Gaussian distribution, while the RMT ap- 
proach can be used in deriving the FGR decay. On the 
other hand, for t relatively long (but shorter than ts): 
the semiclassical approach works well, while the RMT ap- 
proach may meet the problem of deviation of the LDOS 
from the Lorentzian form in the tail region, which is a 
result of the finite domain of the quasi-energy spectrum 

When the underlying classical dynamics has weak 
chaos, non-FGR decay of fidelity may appear in the ex- 
pected FGR regime, due to obvious deviation of the dis- 
tribution of AS from the expected Gaussian distribution 
psf. For this kind of systems, the RMT approach does 
not give a correct prediction for fidelity decay, while the 
semiclassical approach still works. 

In the semiclassical approach, one can separate the 
average fidelity into a mean-value part and a fluctu- 
ating part, denoted by Ma{t) and Mf{t), respectively, 
M(t) = |m(i)p ^Mait) + Mf{t), where 



with average performed over initial states. In the FGR 
regime, the average fidelity is approximately given by the 
mean- value part Ma{t), with M/(i) < Ma{t) 28]. The 
mean- value part Ma{t) can be expressed in terms of the 
distribution P{AS) of the action difference AS, 



Mait) 



where 



(53) 



P{AS)^ ^ / rfrorfpoMAS-AS(po,ro;t)](54) 

In case of weak chaos, a general analytical expression 
for the distribution of AS is still absent. Since Gaussian 
distribution is invalid in this case, it is natural to study 
the possibility of Levy distribution. Due to its infinite 
variance. Levy distribution cannot describe the distribu- 
tion of AS in the long tail region. Therefore, we focus on 
the central part and short tail region of the distribution 
of action difference, which give the main contribution to 
the mean- value part Ma{t) of fidelity. 

We consider the following asymmetric form of Levy 
distribution IItII. 



1 /"^ 

L{x,a,(3) = — / exp{izx)^p{z)dz, 
27r J„oo 



(55) 



Ma{t) = \m{t)\' 



(52) 



with X — AS/e. Here the function V'(^) is 
7/1(2) = exp{—igz — Z)/|z|"[l + i/3 sgn(z) uj{z, a)]}, (56) 
where Lu{z,a) = tan(7ra/2) for a ^ 1, (57) 

tj(z,a) = (2/7r)ln|z| for a =1.(58) 

The parameter a, with < a < 2, determines the decay 
of long tails, i.e., L{x) ~ |a;|~'^^"'""^ for large \x\. The pa- 
rameter P has the domain [—1, 1], with f3 = giving the 
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FIG. 12: Distribution P[(AS - {A5'))/e] of the classical ac- 
tion difference AS at f = 10, for K — 0.4 in the sawtooth 
map, calculated by taking randomly 10^ initial points in the 
phase space. Here {AS) = et{V) = —Tv^et/G, with average 
performed over the phase space. The solid curve is a fit given 
by the Levy distribution in Eq. II55L with a = 1, and /3 and 
Di as two fitting parameters. 



VI. THE LYAPUNOV REGIME 

Increasing the perturbation strength further, one en- 
ters into the Lyapunov regime, in which the average fi- 
delity has the Lyapunov decay e""^* in the special sit- 
uation with negligible fluctuation in the finite-time Lya- 
punov exponent. In the general situation, the fluctuation 
of the finite-time Lyapunov exponent is not negligible, 
and the average fidelity has a Ai(t) decay [s^l, which is 
usually different from the Lyapunov decay. The Ai(t) de- 
cay will be discussed briefly in Sect. IVL^ starting from 
the modified semiclassical approximation to fidelity in 
Eq. 123) • 

As mentioned in the introduction, in the kicked top 
model that has large fluctuation in the finite-time Lya- 
punov exponent, numerical results show that the average 
of the logarithm of fidelity has roughly the Lyapunov de- 
13. In Sect. 



cay [iy. m sect. rVTBl we explain this phenomenon by 
using the technique developed in 30]. 



Ai(t) decay of average fidelity in the deep 
Lyapunov regime 



symmetric distribution, g gives a shift along the x direc- 
tion, and Di is related to the width of the distribution. 
If the Levy distribution can be used as an approxima- 
tion for P{AS), substituting Eq. (jSH) into Eq. one 
obtains. 



Mait) cx exp(-2A ^"), 



(59) 



with the time dependence given by that of Di . Note that 
a Gaussian form of P{AS) corresponds to a = 2, giving 
the well known dependence on a in the PGR decay. 

It is known that the sawtooth map has weak chaos at 
K < 1 and has Cantori structures at K small IsB) . Non- 
FGR decay has been observed in this model, which can 
still be described by Eq. [28| . Therefore, we use this 
model to study the possibility of using the Levy distri- 
bution as an approximation to the distribution P(AS). 

Numerically we use —-ft to fit InM(i), in order to cal- 
culate the decay rate 7 of fidelity. Variation of 7 versus a 
is presented in Fig. for some values of K between 0.1 
and 0.9, before the Lyapunov regime is reached. (The fi- 
delity has been found to have the Lyapunov decay in the 
Lyapunov regime, for K < 1 [l^.) As seen in the Figure, 
for K = 0.1,0.2, and 0.4, there exist some regions of a, 
respectively, in which 7 increases approximately linearly 
with (7. According to Eq. (|59|l . this implies that a ~ 1, if 
the Levy distribution can be used as an approximation to 
P{AS). Figure [T^ shows a fit of the Levy distribution to 
the central part and short tails of the distribution P{AS) 
at K = 0.4, with a = 1 fixed and Di and /? used as two 
adjusting parameters. The agreement is encouraging, for 
which an analytical explanation is still not yet available 
and deserves a further investigation. 



Since Eq. (|26|l can be obtained from Eq. H20|l by replac- 
ing h/£^ with Wp = hD/^, generalization of the results in 
|30l| is straightforward. In this Section, we present the 
main points in the generalization, because part of them 
wiU be used in Sect. lyTBl 

For systems in which AS" has stationary points with 
kp = 0, we denote by a the stationary points of AS and 
by poa the momenta at which kp = 0. When cr 3> 1, sta- 
tionary phase approximation can be used in calculating 
the right hand side of Eq. H2t)|l , which gives 



ruscit) 
mait) 



^mQ(i), where (60) 

a 

V2ihexp [j-AS{poa,r„;t) - (poa -po)Vwp] 



with A^^ 



d^AS{po,ro;t) 



(61) 



Note that Wp in Eq. (|61|) takes the value at the stationary 
point a, with D being a function of poa- 

Let us consider time t > T2, for which there are one 
or more stationary points within the effective window in 
the po space for the integration on the right hand side of 
Eq. (|^ . Similar results can also be obtained for t < T2, 
when the stationary phase approximation is applicable, 
as discussed in Is^. The average fidelity in the ordinary 
sense, with average performed over both Fq andpo, can be 
calculated by using diagonal approximation, when cr 3> 1 
[28l |. and the following result is obtained 



drt) 



Vs 



(62) 
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where Vs ■= Uq-^q- Here, Aa denotes the region 
[Poa'Poa-6] \J[poa+6,Poa]i wheie p^^ = (pOa+PO,Q-l)/2, 
Poa ~ (Poa +Po,a+i)/2, with S being a smaU quantity. 

The main contribution to the integral on the right hand 
side of Eq. H62|l conies from small values of \kp\ in the 
region Vs- For po e Vs close to a stationary point poa, 
kp in Eq. (|29|) can be approximated by 



dt' 



Qj.12 



dr' 
dpo 



dV d^r' 
dr' dpi 



{Po-Paa)- (63) 



Due to exponential divergence of neighboring trajecto- 
ries in phase space, the main contribution to the right 
hand side of Eq. H63() comes from times t' w t. The time 
evolution of the quantity inside the bracket in Eq. I|63|) 
is given by the dynamics of the system described by Hq. 
On average it increases as [Sx{t) / Sx{0)]'^ , with 6x{t) being 
distance in phase space. With increasing time, the num- 
ber of stationary points increases exponentially, roughly 
in the same way as Sx{t)/6x{Q), since the oscillation of 
AS* is mainly induced by local instability of trajectories. 
Then, substituting Eq. into Eq. JHH), we have 



Mit) oc 



1 



_D\Sx{t)/Sx{0)\_ 
When D changes slowly with po and rp, we have 
M{t) cx lA{t) = e-^i(*)*, with 



(64) 



Ai(t) 



lim — In 

Sx(0)^0 t 



5x{t) 



5x{Q) 



(65) 



In systems with constant local Lyapunov exponents, 
Eq. (|65|l reduces to the usual Lyapunov decay with 
Ai(t) = A. On the other hand, when fluctuation in local 
Lyapunov exponent cannot be neglected, coincides 
with the e~'*'i* decay in Ref. f2^ in the limit t oo, with 
Ai = limf^oo Ai(i). 



B. Decay of the average of InM(t) 

To understand the decaying behavior of the aver- 
age of h\M{t) for initial Gaussian wave packets, we 
divide the time t into four time intervals, specifically, 
(0, Ti), (ri,T2), [t2, td), and [td, tg). Here td is a time scale 
defined below, beyond which diagonal approximation can 
be used before average is performed, and tg is the time 
at which the saturation value of fidelity is reached. 

Within the first time interval, M{t) is described by 
Eq. H35|) for most of initial states. For a not large, \nM{t) 
is close to zero. On the other hand, for quite large cr, the 
average of In M{t) can be quite small (cf. Fig.^J. 

For the second time interval (ti,T2), Eq. H39|l can be 
used to calculate InM(t), for initial states with kp satis- 
fying \(Tkp\ S> T^/wp. For this part of initial states, \kp\ 
typically increases as 5x{t) / 5x{Q) [see Eq. (f^ ]. as a re- 
sult, the average of In M(t) for this part of initial states 




FIG. 13: Variation of A(t) and Ai(t) with t in the kicked 
rotator model at A' = 7, with average performed over 10^ 
random initial positions in phase space. The definitions of 
A(t) and Ai(t) are given in Eqs. (1311 and 1651 . respectively. 
The value of A(t) approaches the Lyapunov exponent A ~ 1.27 
quickly, as t increases. 



behaves as — 2A(i)i, where A(t) is defined by Eq. (jJIJ. On 
the other hand, for initial states with small \kp\, which 
lie in the neighborhoods of the stationary points of AS", 
fidelity decay is slower. Since the total size of the regions 
of Po with small \kp\ is small compared with the domain 
of poj the average of In M{t) over all possible initial states 
has a decay rate between h.{t) and 2A{t). In a classical 
system with strong chaos, A(i) approaches the Lyapunov 
exponent A quickly, as seen in Fig. 1131 therefore, in this 
time interval, the decay rate of the average of In M (t) is 
usually close to, or a little larger than A. 

For t around T2 or longer, the main contribution to 
the average of InAf(t) is given by the neighborhood of 
stationary points and one can start from Eq. (|60|l in cal- 
culating M{t). In the third time interval [T2,id), the 
number of stationary points within the effective window 
in Po is small, hence, diagonal approximation cannot be 
used in calculating the absolute value square of the right 
hand side of Eq. H6U|) . We do not know much about the 
decay of the average of In Af (i) in this time interval. One 
should note that the third time interval is quite short, as 
a result of the exponential increment of the number of 
oscillations of AS*. 



In the fourth time interval [trf,^^), diagonal approxi- 
mation is applicable to the absolute value square of the 
right hand side of Eq. H6U|I before average is performed, 
due to the large number of stationary points within the 
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effective window in pq. Then, Eqs. (|60|l and H61|) give 

a. 

2h exp [-2{poa -_po)^/Wp] 



E 



w„ 



(66) 



Using Eq. we write AS"^' as 

^d^V /dr'V dVd^r' 



AS'' 



e I dt' 
'o 



dr'^ \dpt 



dr' dpQ 



(67) 



Arguments similar to those leading from Eq. H63|l to 
Eq. H64|l show that, when the fluctuation of \dx{t) /Sx{0)\ 
is small for po within the effective window in the po space, 
the main decaying behavior of a single fidelity is typically 



M{t) cx \Sx{t)/Sx{ 



(68) 



with Sx(0) being a small displacement from (ro,po)j the 
center of the initial Gaussian packet. Then, it is ready 
to obtain 



exp[ InM(i) ] cx e 



-A(t)t 



-At 



(69) 



where the second equation is obtained, since A(t) is usu- 
ally close to the Lyapunov exponent A within this time 
interval. 

Since the fluctuation of \5x{t) / 5x{Q)\ within the effec- 
tive window in po increases with time, Eq. H69I) becomes 
invalid for t long enough. In the long time limit, the fluc- 
tuation of \5x{t) / 5x{0)\ within the effective window in po 
has similar properties as in the whole po domain, then, 
the average of InM(i) has the Ai(i) decay in Eq. (|65|l . 
with Ai(i) taking the value Ai. For intermediate times, it 
is reasonable to expect that the decaying rate of the av- 
erage of In M(t) decreases from the Lyapunov exponent 
and approaches Ai(t) with increasing t. 

Combining the above results, it is seen that there in- 
deed exists a certain short time interval, in which the 
average of InM(t) follows roughly the Lyapunov decay, 
as observed in the kicked top model [13 ■ Specifically, 
the decaying rate of the average of h\M{t) is close to or 
a little larger than the Lyapunov exponent in the second 
time interval; then, beyond the short third time interval, 
it decreases from the Lyapunov exponent and approaches 
Ai(t) at long time. 

We have tested these predictions, as well as those in the 
previous section, in the kicked rotator model a,t K = 7 , 
which are shown in Fig. 1141 It is seen in the figure that 
the average of InM(t) decays a little faster than the A(t) 
decay (as well as the Lyapunov decay) initially, then, af- 
ter a transient time, it decays a little slower than the 
Lyapunov decay, but obviously faster than the Ai {t) de- 
cay. The predicted long time Ai {t) decay for the average 
of InM(i) is not seen at if = 7, because K is not large 
enough for the value of N taken. We mention that the Ai 
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FIG. 14: Fidelity decay in the kicked rotator model &t K = 7, 
N = 2^"^ , a — 20. Average is performed over 500 initial Gaus- 
sian wave packets with ^ — whose centers are chosen 
randomly in the chaotic sea in phase space. The A{t) and 
Ai{t) decays are the predictions in Eqs. H69|l and H65^ . re- 
spectively, with the values of A{t) and Ai {t) shown in Fig. 1131 
Decay of the average of M{t) is described by the Ai{t) decay. 
The average of \nM{t) is a little faster than the A{t) decay 
in an initial short time interval, as predicted by the theory. 



decay at long time has indeed been observed for K — 10 
in Ref. 20J. 

Finally, for quite large a, the initial approximate Lya- 
punov decay discussed above may disappear. Indeed, in 
this case, | lnM(t)| can be quite large for t < ri, as seen 
in Eq. (|35|l . The average of \nM{t) for short time t is in 
fact dominated by contributions described by Eq. ij^S)), 
even beyond r 1 , since the value of n fluctuates with re- 
spect to the initial condition, i.e., for some initial states ri 
can be relatively long. Hence, the average of In M (t) can 
be obviously smaller than the prediction of the Lyapunov 
decay in the second time interval, as shown numerically 
in 



VII. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have improved the uniform semiclas- 
sical approximation to fidelity by considering the second 
order term in the Taylor expansion of action, which is im- 
portant for fidelity of initial Gaussian wave packets with 
width of the order Vh. Short-time-decay of fidelity is 
analyzed, which is initial-state-dependent, in particular, 
two time scales have been introduced and studied in de- 
tail for initial Gaussian wave packets. Initial FGR decay 
of fidelity in the perturbative regime is confirmed by di- 
rect numerical calculations. Non-FGR decay in the FGR 
regime in systems with weak chaos in the classical limit is 
explained by relating the distribution of action difference 
to the Levy distribution. The average of the logarithm 
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FIG. 15; Schematic diagram for fidelity decay of initial point 
sources in systems with constant local Lyapunov exponents in 
the classical limit, ts is the breakdown time of semiclassical 
approximation, and tH is the Heisenberg time, beyond which 
the Gaussian decay sets in. ad and ar are two perturbation 
scales introduced in Ref. |2^ . Below ad, one has the FGR 
decay for t < ts, and above ar is the Lyapunov regime, in 
which the Lyapunov decay appears beyond the first kick. 

of fidelity is shown to have an approximate Lyapunov de- 
cay within some time intervals, in the Lyapunov regime 
in systems possessing large fluctuations of the finite-time 
Lyapunov exponents in the classical limit. 

As we have demonstrated that fidelity has a decay- 
ing behavior richer than the simple picture mentioned in 
the beginning of the introduction with just four distinct 
regimes. In Fig. 1151 we present a schematic diagram for 
the present understanding of fidelity decay of initial point 
sources in systems possessing constant local Lyapunov 
exponent in the classical limit. For initial Gaussian wave 
packets, short time decay is more complex than for ini- 
tial point sources, for which two time scales ri and T2 
should be introduced, as discussed in this paper. When 
the underlying classical dynamics has large fluctuation 
in the finite-time Lyapunov exponent, the decaying rate 
of fidelity in the Lyapunov regime is not given by the 
Lyapunov exponent. 

Acknowledgments 

The work was supported in part by a Faculty Research 
Grant of National University of Singapore, the Temasek 
Young Investigator Award of DSTA Singapore under 



Project Agreement POD0410553 (BL), and the Natural 
Science Foundation of China No. 10275011 (WGW). 
APPENDIX A: DEPENDENCE OF FIDELITY 
DECAY ON PERTURBATION STRENGTH IN 
DEEP LYAPUNOV REGIME FOR INITIAL 
POINT SOURCES 

In this appendix, we consider the dependence of fidelity 
decay on perturbation strength in the deep Lyapunov 
regime for initial point sources, which can be estimated 
by using Eq. (|21|l . Let us first divide the domain [0, 2ti) of 
Po into a series of subregions, [poj,Po(i-i-i))) so that within 
each of the subregions the phase AS/h of the integrand 
on the right hand side of Eq. (PT|l can be approximated 
by a linear function, 

AS/n ~ a{kjPo + bj), for po e boj,Po(i+i)), (Al) 

where the parameters kj and bj do not depend on pq. 
Here pqj = for j — 1 and poj = 2ti for j — Nx + 1, with 
Nx denoting the number of the subregions. Substituting 
Eq. lPT|l into Eq. (EIJ, we have 

^ Nx 

mn{t) ~ > Xj with 

i=i 

Xj ■= — (^g'^'^'^iPoU + i) _ ^i'^^jPOj^ ^icrbj (A2) 

kj 

In arranging the subregions, we require that Nx should 
be as small as possible, conditional on that the linear ap- 
proximation in Eq. I|A1|I does not lose the main contribu- 
tion to the right hand side of Eq. (|21|l . Since the phase 
AS/h is proportional cr, Nx increases approximately lin- 
early with a. In the deep Lyapunov regime with quite 
large a, for the subregions chosen in this way, the phase 
of Xj, which is approximately proportional to cr, can be 
regarded as random with respect to j and tq, then, diago- 
nal approximation can be used in computing the averaged 
fidelity, 

Nx 

M,[t) = ^1 dro|m,(OP ^J-^^^, (A3) 

where we assume that fluctuation of Nx with respect to 
To is small. Finally, using Nx oc cr, we have the following 
approximation, 

Mp(i) cx 1/ct, for cr > 1. (A4) 
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